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Abstract 



Cosmological models with two interacting fluids, each satisfying the 
strong energy condition, are studied in the framework of classical Gen- 
eral Relativity. If the interactions arc phcnomenologically described by 
a power law in the scale factor, the two initial interacting fluids can be 
equivalently substituted by two non interacting effective fluids, where one 
of them may violate the strong energy condition and/or have negative 
energy density. Analytical solutions of the Friedmann equations of this 
general setting are obtained and studied. One may have, depending on 
the scale where the interaction becomes important, non singular universes 
with early accelerated phase, or singular models with transition from de- 
celerated to accelerated expansion at large scales. Among the first, there 
are bouncing models where contraction is stopped by the interaction. In 
the second case, one obtains dark energy expansion rates without dark 
energy, like ACDM or phantomic accelerated expansions without cosmo- 
logical constant or phantoms, respectively. 

PACS number(s): 98.80.Cq, 98.80.-k, 98.80.Bp 
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1 Introduction 



In general, the fluids describing the matter content of a cosmological model are 
considered to be non interacting. For instance, in the Standard Cosmological 
Model, cold dark matter, neutrinos, dark energy and the baryon-photon fluid 
are usually taken to be decoupled. However, it is well known that in the early 
Universe many of these fluids were coupled through annihilation and/or scat- 
tering processes [1]. Many approachs to describe these primordial interactions 
have been proposed [2, 3, 4, 5, 6, 7]. Furthermore, interactions between the 
two dark components have been investigated in order to explain the coinci- 
dence problem, and effective phantom acceleration without phantomic equation 
of state [8, 9, 10, 11]. Finally, interactions between two dark matter compo- 
nents, or dark matter self-interaction [12], and interaction between dark matter 
and baryons [13] have also been investigated. In the last case, the possibility of 
interactions between baryons and dark matter due to the strong force has been 
proposed because it can solve some conflicts between numerical simulations on 
structure formation with observations. One can justify the fact that such strong 
interactions have not yet been detected because most of the WIMP detectors 
are ground based [13]. 

In this report, we explore the consequences of some phenomenological mod- 
els of interaction between two positive energy fluids with constant equation of 
state, both satisfying the strong energy condition. We follow the line of research 
of Refs [10, 11], postponing the physical justification of these phenomenological 
proposals for future publications in order to concentrate on their physical conse- 
quences, and their potentialities to solve cosmological puzzles (e.g, concordance 
problem, cosmological singularities), with the hope that they may be tested soon 
[14] . Assuming a power law form in the scale factor for these interactions, similar 
to what was done in Ref. [11] but not restricting to dark energy-dark matter in- 
teractions, we arrive at an effective Friedmann equation containing two effective 
non interacting perfect fluids where one of them, depending on the interaction, 
may have arbitrary equation of state, including phantom-like behaviour, and/or 
negative energy density^. This gives rise to many types of cosmological analyt- 
ical solutions which we obtain explicitely: bouncing universcis, as those studied 
in Refs. [16, 17, 18], and singular models with late acceleration phases, like the 
ACDM model or others with phantomic acceleration leading to a big-rip. In 
all these cases, the occurrence of interactions between the fluids may lead, even 
in the framework of classical General Relativity, to accelerating phases in the 
Universe, without the need to suppose that such fluids violate the strong energy 
condition. 

The paper is organized as follows: in Section II we describe our phenomeno- 
logical models and we arrive at the effective Friedmann equations, from which 
we obtain the analytical solutions presented and discussed in Section III. Section 

IV is devoted to the conclusions. 

^Some similar situations were studied in Refs. [15]. 
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2 Cosmology with two interacting fluids 

We model the interaction between two fluids satifying the strong energy condi- 
tion and with equations of state Pi = WiPt, i € {1,2}, with Wi =const.> —1/3, 
and Pi > as follows: 

Pi+3H{l + wi)pi = -Q, 

P2+3H{1+W2)P2 = Q, (1) 

where H = a/a \s the Hubble parameter, Q is the interaction rate, and a dot 
represents derivative in cosmic time. 

We assume, without loss of generality, that the energy density of the second 
fluid can be written as: 

/ /ao\3(i+«'2) 
P2=P2o-^( — ) , (2) 

JO ^ ^ 

where p2o , /o and oq are constants, and / is an arbitrary time-dependent func- 
tion. It follows from (1) and (2) that 



Q ^ P2- = P2o— [—) (3) 



We now assign an ansatz for /, 

/o V«o 



/ / a ^ -'"^ 



(4) 



where w/ is a constant. Hence 

[-3(1 + 102+^/)] 

Consequently 



p2=p2o[^\ ■ (5) 



/ \ -3(l+W2+l«/) 

Q = -3wfP2oH i^-j , (6) 

The solution of the conservation equation for the first fiuid 

/ ^ \ -3{1+W2+Wf) 

pi+3H{l + wi)pi=3wfP2oH l—j (7) 

is given by 

(\ —3(1+W2+Wt) 
-) , (8) 

do J 

with Ci a constant. The total energy density then reads 

/ \ / \ — 3(l+J«2+«'f ) 

p^ = Cia-^(i+-) + ^^"^"'^'"^^ ( ^\ . (9) 

w\—W2— Wf \ao J 
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Making the definitions, 

■:Ul+<r2 + u-.) , X 

W+=Wi, W-=W2+Wf, p+ = Ci, p- = '-^^ '-^ ^ 10 

W+ — W- 

the Friedmann equation can be written as: 

= (p+a-3(i+-+) + p_a-3(i+«'-)) , (11) 

where = SttG/S. 

Hence, from two interacting fluids satisfying the strong energy condition, 
Wi,W2 > —1/3, and positive energy density, pi,p2 > 0, interacting via Eq. (1) 
with Q given in (6), one obtains two effective non interacting fluids characterized 
by the parameters w+,W-, p+, p- in the Friedmann equation (11), where one 
satisfies the strong energy condition, Wi = w+ > —1/3, but the other can 
violate the strong energy condition and/or have negative energy as long as, 
from deflnitions (10), and depending on the values of Wf in (6), one may obtain 
W- < —1/3 and/or p- < 0. This fact leads to interesting cosmological solutions 
where bounces and/or late accelerated phases may happen. We present these 
possibilities in the following section. Note that if the two fluids have the same 
equation of state, Wi = W2 there is no second effective fluid. Also, if w+ = 
the quantity p_ cannot be defined. Hence we do not consider this singular point 
in parameter space below. 



3 Analytical solutions 

By introducing a new coordinate time r: 



dT=^, with P =^{2w+ -w- + 1) (12) 

CI ^ 



we have as solution for the scale factor: 



with 



a(r) = Ob ( — - — ) , (13) 
^0 P- 



a = 777 T (14) 

3(u;_ -w+) ^ ' 

With the new coordinate time r, it is therefore possible to generalize the 
solution in the conformal time found in [17], and the cosmological models with 
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dust plus dark energy with arbitrary constant equation of state obtained in [20] , 
to arbitrary values of 'W-^- , W- . 
The Hubble function is 

a 2aT 



while cosmic acceleration reads 
d 2a 
a 



(1 + 3«;+)q^ + ^ 
■^0 P- . 



(17) 



(18) 



One may have transitions from decelerated to accelerated phases, and vice- versa, 
when 

_2 \P-\ T-,7 



, - , ■. (19) 

p- a{l + 3w+) ^ ' 

The asymptotic behaviours occur for r — > ±oo and, if p_ > 0, for r ±to. 
In the first case one has 

t oc |r|3«(i+»+), a(t) oc |t|2/[3(i+«'+)]^ (20) 
while for the second one obtains, 

i oc |t ± ro|3«(i+--)/2, a{t) oc |f|2/[3(i+--)], (21) 

Finally, looking at the evolution of the interaction term (6), one should 

expect that its modulus \Q\ would decrease when the Universe expands and 
Hubble time tn = ^/H increases, which is indeed the case of Eq. (6) unless 
w_ < —1 ^ Wf < —1 — W2- However, this extreme situation can be explained 
using Eqs. (5,8), where it is shown that the energy densities of the two fluids pi 
and p2 increase with expansion in that case, making the interaction stronger. 
We will analyze these situations below in more details using the equation 



2a 



0/2 I. I ^2(l+«;_)/(«;_-«;+) 



fr^ \p-\y 



3a(3 + 2w- + w+)^ + — 
^0 P- . 



.(22) 



We will now discuss in more details the possible cases in terms of the sign 
of As we are assuming that both original fluids satisiy the strong energy 
condition, then w+ = wi > —1/3. 

3.1 p_ < 

The scale factor reads: ^ ^ 

a(r) = ab(^^ + l^ , (23) 

implying that —00 < r < 00, and hence, from Eq. (17), the scale factor has an 
extremum at r = 0. 
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3.1.1 a>0 



When a > 0, Eq. (18) shows that this cxtrcmum is a minimum, and wc have 
a non singular bouncing universe. There are transitions from deceleration to 
acceleration and vice-versa because there are real roots from Eq. (19) in that 
case: a(l + 3w+) > 0. As w_ > w+, the negative energy fluid dominates when 
the universe is small, avoiding the singularity in the accelerated phase, while 
the positive energy fluid dominates when the universe is large, which expands 
deceleratcly in this regime, as it can be seen from Eq. (20). 

Regarding the interaction term Q, one can see from Eq. (22) that for 

l! 1 

°< t2 < 3a(3 + 2«;_+u;+) 

Q is increasing while the model expands. This can be understood by noticing 
that this is the period near the bounce^, where the scale factor increases slowly 
while the Hubble time th = 1/H decreases abruptly, largely compensating the 
slow expansion. 



3.1.2 a<0 

When a < 0, Eq. (18) shows that the cxtrcmum is a maximiun. When the 
universe is small, the w+ fluid dominates and the W- fluid becomes important 
only around the maximum. We have a big-bang big-crunch model (remember 
wc are assuming > —1/3, hence t ±oo implies t finite in this case, see 
Eq. (20)). There are no transitions from deceleration to acceleration and vice- 
versa because a(l -|- 3w+) < in this case (see Eq. (19)), whatever is the value 
of w_ . 

The interaction term may increase with expansion only for the period shown 
in Eq. (24), now around the maximum of the scale factor, if W- < — (3-I-W+) /2 < 
—4/3 (remember that w_ < w+ in this case). This can be justified in terms of 
the slow increasing of a in that period, the increasing of the energy densities of 
the two fluids pi and p2 because W- < — 1. 



3.2 p_ > 

In this case the scale factor reads 

air) = Ub (^^ - ?j , (25) 

implying that — oo < r < — tq, or tq < t < oo. Hence, from Eq. (17), the scale 
factor has no extremum being either an always contracting or always expanding 
model. We will concentrate on the always expanding solutions. 

The conditions for having transitions from deceleration to acceleration and 
vice-versa in this case are (see Eq. (19)), besides a{l + 3w+) < 0, the condition 

*As 3 + 2w- + w+ > 2 in this case, [3a(3 + 2w- + w+)]-^ in Eq. (22) can be big only if 
< q; << 1, but then the bounce also lasts very long in r (see Eq. (19)) 
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that \t\ > \to\ (sec Eq. (13)), which implies that |a(l + 3w+)| < 1. Hence, in 
this case, a > imposes that w+ < — 1/3 and W- > —1/3, and a < imphes 
that w+ > — 1/3 and W- < —1/3, as it should be if the two effective fluids have 
positive energy. As we are assuming w+ > —1/3, there are no transitions when 
a > 0. 

3.2.1 a > 

When a > 0, both effective fluids satisfy the energy conditions and have positive 
energy. These are the standard cases of two ordinary fluids governing a universe 
with decelerated expansion from an initial singularity. The interaction term Q 
always decreases with expansion. 

3.2.2 a < 

When a < 0, the W-^- fluid dominates when the universe is small, necessarily 
reaching a singularity (from Eq. (20), r — > ±oo implies t finite in this case), and 
the w- fluid dominates when it is large. If w_ < — 1/3 one can have a transition 
from decelerated to accelerated expansion, and if W- < — 1, one gets a big rip. 

As pointed above, in this last case one may have increasing \Q\ with ex- 
pansion. Looking at Eq. (22), one can see that, for w_ < — (3 + w+)/2, or 
Wf < —{3 + Wi+ 2w2)/'2, one has the bizarre situation where \Q\ always in- 
creases with expansion, even near the initial singularity, when the ordinary 
= wi fluid dominates and expansion takes place in the standard way. It 
seems that the increasing in p2 is so strong (see Eq. (5)), that it compensates 
expansion, and the decreasing of p- and H. When —(3 + w+)/2 < W- < — 1 
one has the more reasonable situation where \Q\ decreases with expansion when 
the ordinary w+ = wi fluid dominates and increases only in the period given by 



when we expect that the energy densities of the two fluids pi and p2 increase 
with expansion (see Eqs. (5,8)), making the interaction stronger, and Hubble 
time tn = ^/H decreases. Note that, depending on the value of and W-, 
the increasing of |<5| niay happen after p-p+ equilibrium, which takes place at 

= 2tq. For instance, if w+ = 0, there is a period when \Q\ decreases and the 
W- fluid dominates if —6/5 < w;_ < — 1. 

For w- > —1, IQI always decreases with expansion. 

3.3 Final remarks 

For the sake of completeness, let us mention some different models one may 
obtain if the condition w+ > —1/3 is relaxed. One may have pre-big bang [19] 
and singular inflationary models without graceful exit when p_ > 0. When 
P- < 0, there are bouncing models between big rips or between accelerated 
contraction and expansion, and models which expands from a pre-big bang like 




3a{5 + 2w-+w+)' 



1 



(26) 
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initial state, till a maximum size, and then contracts to a time reversed pre-big 
bang like final state, with transitions from accelerated to decelerated phases and 
vice-versa around the maximum. These are not realistic models. 

From all these possibilities, the most interesting cases are the non singular 
solutions and the models with late accelerated expansion. Among the first, 
there are the bouncing models with p_ < and a > 0. They are non singular 
models which are connected to a standard expansion phase dominated by a fluid 
satisfying the strong energy condition when the universe is large. In particular, 
the bouncing model studied in Ref. [16] can be obtained with the choice wi = 
1/3, W2 = and wj = 1, yielding an interaction term given hy Q = —3KHa~^, 
which is strongly suppressed when a is large. One can view this interaction 
as happening in a temperature where baryons are relativistic but dark matter 
is not, and they interact through the strong force as suggested in Ref. [13]. 
This could also happen in a temperature where both are relativistic, but not 
exactly with the same equation of state parameter. In that case on should have 

« 2/3 and g « ^2KHa-^. 

In the second group there are the models with transition from decelerated 
to accelerated expansion, as in the cases with parameters p_ > 0, a < 0, 
and 7i>_ < —1/3. One possibility would be to have two different dark matter 
components, or dark matter and baryons, with 0<u'i<<l,0<W2<<l and 
wi > W2, interacting with Wf = —1, yielding an effective ACDM model with 
w;+ w and « — 1. In this case, one has Q « 3KH, a mild decreasing with 
expansion depending only on the Hubble parameter. Note that, as the effective 
cosmological constant is given by A oc p_ oc {wi — W2), its smallness could then 
be related with the smallness of wi and W2- 

4 Conclusion 

Wc have obtained some interesting cosmological models with flat spatial sections 
from ordinary fluids which interact with the phenomenological interaction rate 
(6). Bouncing models, which cannot be obtained from two non interacting fluids 
within classical General Relativity unless one of them violates the null energy 
condition [21], may arise in this context even if the two interacting fluids can 
be modelled by dust and/or relativistic fluids, like dark matter components and 
baryons, depending on the background temperature. This same combination of 
fluids, with another choice of the power Wf in Eq. (6), may lead to the ACDM 
expansion rate without the need to introduce any fundamental cosmological 
constant, whose small observational value is a challenge for particle physics 
theory. 

Of course the relevance of the phenomenological interaction rate (6) in differ- 
ent phases of the history of the real Universe must be justified in micro-physical 
grounds. We will come back to this problem in forthcoming publications. Also, 
the evolution of primordial quantum perturbations in such models must be car- 
ried out in order to compare the results with CMBR data. For instance, one 
should investigate if the bouncing models presented here could also lead to scale 
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invariant spectra of perturbations as it was verified in other bouncing universes 
[22]. These are subjects beyond the aim of the present report, which is to point 
out that it is maybe not necessary to evoke the existence of exotic and/or un- 
known substances to yield the present accekration of the Universe and/or to 
avoid the initial singularity within classical General Relativity. 
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